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Abstract. Discrete geometric estimators approach geometric quantities 
on digitized shapes without any knowledge of the continuous shape. A 
classical yet difRcult problem is to show that an estimator asymptot- 
ically converges toward the true geometric quantity as the resolution 
increases. We study here the convergence of local estimators based on 
Digital Straight Segment (DSS) recognition. It is closely linked to the 
asymptotic growth of maximal DSS, for which we show bounds both 
about their number and sizes. These results not only give better insights 
about digitized curves but indicate that curvature estimators based on 
local DSS recognition are not likely to converge. We indeed invalidate 
an hypothesis which was essential in the only known convergence the- 
orem of a discrete curvature estimator. The proof involves results from 
arithmetic properties of digital lines, digital convexity, combinatorics, 
continued fractions and random polytopes. 



1 Introduction 

Estimating geometric features of shapes or curves solely on their digitization is 
a classical problem in image analysis and pattern recognition. Some of the geo- 
metric features are global: area, perimeter, moments. Others are local: tangents, 
normals, curvature. Algorithms that performs this task on digitized objects are 
called discrete geometric estimators. An interesting property these estimators 
should have is to converge towards the continuous geometric measure as the 
digitization resolution increases. However, few estimators have been proved to 
be convergent. In all works, shapes are generally supposed to have a smooth 
boundary (at least twice diflerentiable) and either to be convex or to have a 
finite number of inflexion points. The shape perimeter estimation has for in- 
stance been tackled in [TT]. It proved the convergence of a perimeter estimator 
based on curve segmentation by maximal DSS. The speed of convergence of 



several length estimators has also been studied in [4]. Klette and Zunic [TO] 
survey results about the convergence (and the speed of convergence) of several 
global geometric estimators. They show that discrete moments converge toward 
continuous moments. 

As far as we know, there is only one work that deals with the convergence 
of local geometric estimators [3]. The symmetric tangent estimator appears to 
be convergent subject to an hypothesis on the growth of DSS as the resolution 
increases (see Hypothesis |4T|) . The same hypothesis entails that a curvature 
estimator is convergent: it is based on DSS recognition and circumscribed circle 
computation (see Definition [9]) . 

In this paper, we relate the number and the lengths of DSS to the number 
and lengths of edges of convex hulls of digitized shapes. Using arguments related 
to digital convex polygons and a theorem induced by random polytopes theory 
[l], we estimate the asymptotic behaviour of both quantities. We theoretically 
show that maximal DSS do not follow the hypothesis used in |3j. Experiments 
confirm our result. The convergence theorem is thus not applicable to digital 
curves. As a consequence, the existence of convergent digital curvature estimators 
remains an open problem. The paper is organized as follows. First, we recall some 
standard notions of digital geometry and combinatoric representation of digital 
lines, i.e. patterns. The relations between maximal segments and edges of convex 
digital polygons are then studied to get bounds on maximal segments lengths and 
number. Finally, the asymptotic behaviour of maximal segments is deduced from 
the asymptotic behaviour of convex digital polygons. The growth of some DSS 
is thus proved to be too slow to ensure the convergence of curvature estimation. 
This theoretical result is further confirmed by experiments. 

2 Maximal digital straight segments 

We restrict our study to the geometry of 4-connected digital curves. A digital 
object is a set of pixels and its boundary in is a collection of vertices and 
edges. The boundary forms a 4-connected curve in the sense used in the present 
paper. Our work may easily be adapted to 8-connected curves. In the paper, all 
the reasoning are made in the first octant, but extends naturally to the whole 
digital plane. The digital curve is denoted by C. Its points [Ck) are assumed to 
be indexed. A set of successive points of C ordered increasingly from index i to 
j will be conveniently denoted by [CiCj] when no ambiguities are raised. 

2.1 Standard line, digital straight segment, maximal segments 

Definition 1. (Reveilles \14V The set of points {x,y) of the digital plane veri- 
fying ^ < ax ~ by < ^ + \a\ + \b\, with a, b and ji integer numbers, is called the 
standard line with slope a/b and shift /i. 

The standard lines are the 4-connected discrete lines. The quantity ax — by is 
called the remainder of the line. The points whose remainder is /i (resp. + 



\a\ + \b\ — 1) are called upper (resp. lower) leaning points. The principal upper 
and lower leaning points are defined as those with extremal x values. Finite 
connected portions of digital lines define digital straight segment. Since we work 
with restricted parts of C, we always suppose that indices are totally ordered on 
this part. 

Definition 2. A set of successive points [CiCj] of C is a digital straight segment 
(DSS) ijj there exists a standard line D{a,b,fi) containing them. The predicate 
"[CiCj] is a DSS" is denoted by S{i,j). 

The first index j, i < j, such that S{i,j) and -^S{i,j + 1) is called the front 
of i. The map associating any i to its front is denoted by F. Symmetrically, the 
first index i such that S{i,j) and ^S{i — is called the back of j and the 
corresponding mapping is denoted by B. 

Maximal segments form the longest possible DSS in the curve. They are 
essential when analyzing digital curves: they provide tangent estimations |6|13) . 
they are used for polygonizing the curve into the minimum number of segments 

m- 

Definition 3. Any set of points [C'iCj] is called a maximal segment iff any of 
the following equivalent characterizations holds: (1) S{i,j) and ^S{i, j + 1) and 
-^S{i-l,j), (2) B{j)=i andF{i)=j, (3)3k,i = B{k) and j = F{B{k)) , (4) 
3k',i = B{F{k')) andj = F{k'). 

From characterizations (3) and (4) of Definition O any DSS [QCj] and 
hence any point belongs to at least two maximal segments (possibly identical) 

[Cb(j)Cf(b(j))\ and [CB(F{i))CF(i)\- 

2.2 Patterns and DSS 

We here recall a few properties about patterns composing DSS and their close 
relations with continued fractions. They constitute a powerful tool to describe 
discrete lines with rational slopes |2|8j . Since we are in the first octant, the slopes 
are between and 1. 

Definition 4. Given a standard line {a,b,^), we call pattern of characteristics 
(a, b) the succession of Freeman moves between any two consecutive upper leaning 
points. The Freeman moves defined between any two consecutive lower leaning 
points is the previous word read from back to front and is called the reversed 
pattern. 

A pattern (a, 6) embedded anywhere in the digital plane is obviously a DSS 
(a, b, n) for some /i. Since a DSS contains at least either two upper or two lower 
leaning points, a DSS (a, 6, /i) contains at least one pattern or one reversed pattern 
of characteristics (a, b). 



Definition 5. We call simple continued fraction and we write: 
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We call fc-th convergent the simple continued fraction formed of the k + 1 first 
partial quotients: Zk = ^ = [0,ui, . . . , Uk] . 

There exists a recursive transformation for computing the pattern of a standard 
hne from the simple continued fraction of its slope [2] . We call E the mapping 
from the set of positive rationnal number smaller than one onto Freeman-code's 
words defined as follows. First terms are stated as E{zo) = and E{zi) = O"^! 
and others are expressed recursively: 



In the fohowing, the complexity of a pattern is the depth of its decomposition 
in simple continued fraction. We recah a few more relations: 



We now focus on computing vector relations between leaning points (upper 
and lower) inside a pattern. In the following we consider a DSS (a, 6, 0) in the first 
octant starting at the origin and ending at its second lower leaning point (whose 
coordinate along the x-axis is positive). We define a/b = Zn = [0, ui, . . . , u„] 
for some n. Points will be cahed Ui,Li, U2 and L2 as shown in Fig. [TJ We can 
state UiLi = U2L2 and U1U2 = L1L2 = (6, a). We recall that the Freeman 
moves of [UiLi] are the same as those of [U2L2]. Furthermore Freeman moves 
between Ui and U2 form the pattern (a, 6) and those between Li and L2 form 
the reversed pattern (a, 6). 

Proposition 1. A pattern with an odd complexity (say n = 2i + 1) is such that 

UiLi = (M2i+l-l)(<?2i,P2i) + (92i-l,P2i-l) + (l,-l) anrfLiU2 ((72j-1,P2j + 1)- 

Moreover the DSS [UiLi] has i;(z2i)"''+'"^ « left factor, and the DSS [L1U2] 
has i?(z2i-i)"^' as a right factor. 

Proof. From Eq. ^ we have: P2i+iq2i — P2iq2i+i = (— l)2'+i+i = 1^ which can 
be rewritten as: aq2i — bp2i = 1. {q2i,P2i) are clearly the Bezout coefficients of 
(a, b). One can check that point (6 + 1 — q2i, a — 1 — p2i) is Li: its remainder is 
a+b— 1 and its x-coordinate while positive is smaller than b. We immediately get 
UiLi = (&+ 1 - (?2j, a - 1 -P2i). Using Eq. (01 yields: UiLi = {{u2i+i - l)92i + 
g2i-i + l, (w2i+i-l)p2i+P2i-i-l). FromLiU2 = -U1L1+U1U2, we further get 



E{z2^+l) ^ E{z2^)'''' + 'E{z2^-l) 
E{Z2^) = E[z2^-2)E{z2^-lT'^ 



(1) 

(2) 



Pkqk-i - Pk-iqk = (-1) 

{Pk, qk) = Uk{pk-i,qk-i) + (pfc-2, gfe-2) 



(3) 
(4) 
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Fig. 1. A DSS(a,6, 0) with an odd complexity slope, taken between origin and 
its second lower leaning point. 



that L1U2 = ((72»-l,P2( + l)- FromEq. P £:(z2j)""'+'"^ is a left factor of [;7iC/2] 
but also of [UiLi]. Writing E{z2i+i) as £:(z20""'+'"^£^(^2*-2)£^(z2i-i)""'+\ and 
expanding L1U2 as {u2iq2i~i + 921-2 - 1, U2iP2i-i +P2i-2 + 1) with Eq. we 
see that £'(z2i_i)"^* is a right factor of [Lit/2]. □ 

Proposition 2. A pattern with an even complexity (say n = 2i) is such that 

UlLi = ((72t-l + l,P2i-l-l) andLiU2 = (w2i-l)((72i-l,P2j-l) + (l?2i-2,P2i-2) + 

(— 1, 1). Moreover the DSS [UiLi] has £^(222-2)"^'"^ as a left factor, and the DSS 
[L1U2] has E{z2i-i)^^^~^ as a right factor. 

3 Properties of maximal segments for convex curves 

In this section, we study relations between maximal segments and digital edges 
of convex shape digitization. The dilation of 5 by a real factor r is denoted by 
r ■ S. Let Vm be the digitization of step 1/m, i.e. if 5 is a real shape: VmiS) = 
(m • S") n Z^. The length estimator based on the city-block distance is written as 
£1. 

3.1 Convex digital polygon (CDP) 

Definition 6. A convex digital polygon (CDP) F is a subset of the digital plane 
equal to the digitization of its convex hull, i.e. F = I?i(conv(P)). Its vertices 
iyi)i=i..e form the minimal subset for which F = I>i(conv(Vi, . . . , Vg)). The 
points on the boundary of F form a ^-connected contour. The number of vertices 
(or edges) of F is denoted by ne{F) and its perimeter by Per(J^). 

A CDP is also called a lattice convex polygon [16j. An edge is the Euclidean 
segment joining two consecutive vertices, and a digital edge is the discrete seg- 
ment joining two consecutive vertices. It is clear that we have as many edges as 



digital edges and as vertices. From characterizations of discrete convexity [5], we 
clearly see that: 

Proposition 3. Each digital edge of a CDP is either a pattern or a succession 
of the same pattern whose slope is the one of the edge. In other words, both 
vertices are upper leaning points of the digital edge. 

We now recall one theorem concerning the asymptotic number of vertices 
of CDP that are digitization of continuous shapes. It comes from asymptotic 
properties of random polytopes. 

Theorem 1. (Adapted from Balog, Bdrdny f^) If S is a plane convex body with 
boundary and positive curvature then T>m{S) is a CDP and 

ci{S)m^ < neiV,n{S)) < C2iS)m'i 

where the constants Ci{S) and 02(8) depend on extremal bounds of the curvatures 
along S. Hence for a disc ci and C2 are absolute constants. 

3.2 Links between maximal segments and edges of CDP 

Maximal segments are DSS: between any two upper (resp. lower) leaning points 
lays at least a lower (resp. upper) leaning point. The slope of a maximal segment 
is then defined by two consecutive upper and/or lower leaning points. Digital 
edges are patterns and their vertices are upper leaning points (from Prop. [3]). 
Thus, vertices may be upper leaning points but never lower leaning points of 
maximal segments. We have 

Lemma 1. A maximal segment cannot be strictly contained into a digital edge. 
We now introduce a special class of digital edge. 

Definition 7. We call supporting edge, a digital edge whose two vertices define 
leftmost and rightmost upper leaning points of a maximal segment. 

Relations between maximal DSS and digital edges are given by the following 
lemmas: 

Lemma 2. A supporting edge defines only one maximal segment: it is the only 
one containing the edge and it has the same slope. If a maximal segment contains 
two or more upper leaning points then there is a supporting edge linking its 
leftmost and rightmost upper leaning points with the same slope. If a maximal 
segment contains three or more lower leaning points then it has a supporting 
edge. 

Lemma 3. // a maximal segment is defined by only two consecutive lower lean- 
ing points then it has one upper leaning point which is some vertex of the CDP 
by convexity. 



Lengths of maximal segments and digital edges are tightly intertwined, as 
shown by the two next propositions. 

Proposition 4. Let [VkVk+i] be a supporting edge of slope ^ made of f patterns 
(a, 6) and let MS be the maximal segment associated with it (Lemma\B). Their 
lengths are linked by the inequalities: 

C\VkVk+i) < C'{MS) < :^C\VkVk+i)-2 and ^C\MS) < C'iVkVk+i) < C^iMS) < 3C\VkVk+i) 

Proof. Vertices Vk and Vk+i are leftmost and rightmost upper leaning points of 
MS. The points Vk — (6, a), Vk+i + {b,a) while clearly upper leaning points of 
the standard line going through [VfcT4+i] cannot belong to the CDP. Hence MS 
cannot extend further of its supporting edge of more than \a\ + \b\ — 1 points 
on both sides. Consequently C\MS) < C^{VkVk+i) + 2{\a\ + \b\ - 1). Using 
C\VkVk+i) = /(|a| + |5|) brings: C^{VkVk+i) < C\MS)) < ^^C\VkVk+i) -2. 
Worst cases bring C^VkVk+i) < C^{MS) < iC^{VkVk+i) □ 

Proposition 5. Let MSk' be a maximal segment in the configuration of Lemma\^ 
and so let Vk be the vertex that is its upper leaning point. The length of the max- 
imal segment is upper bounded by: 

C\MSk') <^{C\Vk-iVk)+C\VkVk+i)) 

Proof. We call Li, L2 the leftmost and rightmost lower leaning points and 
U2 = Vk the upper leaning point (see Fig. [J). Suppose that MSk' has a slope 
with an odd complexity (say 2i + l). Proposition [T] implies £^(LiU2) = q2i+P2i- 
There is clearly a right part of [L1C/2] (i-e. [-Z^iV^]) that is contained in [Vfe_iT4] 
and touches Vk. The pattern i?(z2i-i)"^' is a right factor of [iiC/2] (Proposi- 
tion [T] again) . It is indeed a right factor of [V^-iVfe] too, since it cannot ex- 
tends further than Vk-i to the left without defining a longer digital edge. We 
get [Vk-iVk] 3 E{z2^-l)'''-' and immediately C^iVk-^iVk) > U2^C^{E{z2^^l)) = 

U2i{q2i-1 + P2i-l). 

From Eq. (g]), we have: q2i + P2i = U2iiq2i-i + P2i-i) + q2i-2 + P21-2 and 
922-2 +P21-2 < 921-1 +P21-1- We obtain immediately £^(LiU2) = q2i + P2i < 
{u2i + l){q2i-i -\- P2i-i)- By comparing this length to the length of the digital 
edge [Vk-iVkl we get ^^(LiUa) < ^C\Vk-iVk). 

Propositiondland similar arguments on [VfeV^+i] brings £^(U2L2) < ^'^^^^■^ L^ iVk-iVk) . 
Worst cases are then £i(LiU2) < 2C^{Vk-iVk) and £i(U2L2) < 2C^{VkVk+i). 
The case where MSk' has a slope with an even complexity (say 2i) uses Prop. [2] 
and is treated similarly. 

Since MS has only one upper leaning point, it cannot be extended further 
than £^(U2L2) on the left and £^(LiU2) on the right (Lemma[2]). We thus get 
C\MSk') < A{C\Vk-iVk) + C\VkVk+i)). □ 

A proof of the following theorem based on pattern analysis is given in Ap- 
pendix B for limited space reasons. A similar result related to linear integer 
programming is in [T^. It may also be obtained by viewing standard lines as 
intersection of two knapsack polytopes [9] . 



Theorem 2. Let E he a supporting edge whose slope has a complexity n, n >2, 
then the maximal segment containing E includes at most n other edges on each 
side of E. 

Corollary 1. The shortest pattern of a supporting edge for which its maximal 
segment may contain 2n + I digital edge is = [0,2,..., 2]. If the DCP is 
enclosed in a m x m grid, then the maximal number n of digital edges included 
in one maximal segment is upper bounded as: n < log ^/log (1 + -^2) — 1- 

Proof. The number L = [0, 2, . . . , 2, . . .] is a quadratic number equal to — 1 + \/2- 
Its recursive characterization is Un — 2[/„_i + [/„_2 with Uq = and Ui = 1. 
Solving it leads to J7„ = ^ ((1 + \/2)" - (1 - \/2)") . Hence asymptotically, 

C/n « ^(1 + V2)" and lim„^oo ij^.^L. 

The shortest edge of slope complexity n is clearly an n-th convergent of L. 
To fit into an to x m grid, the complexity n is such that Un+i < m. We thus 
obtain that n < log ^/log (1 + ^2) - 1. □ 

Proposition 6. There exists at most two maximal segments per vertices in the 
configuration of Lemma [3 with different parities of complexity. 

Proof. We first prove that there is at most one maximal segment with only one 
upper leaning point on a vertice of a DCP with an even complexity. 

Let us suppose that MSi and MS2 are two maximal segments sharing a 
vertice of the CDF (say U2) with their slopes 



3.3 Asymptotic number and size of maximal segments 

We assume in this section that the digital convex polygon F is enclosed in a 
mx m grid. We wish to compute a lower bound for the number of edges related 
to at least one maximal segment. We show in Theorem [3] that this number is 
significant and increases at least as fast as the number of edges of the DCP 
divided by log to. From this lower bound, we are able to find an upper bound 
for the length of the smallest maximal segment of a DCP (Theorem [4]). We first 
label each vertex of the DCP as follows: (i) a 2-vertex is an upper leaning point 
of a supporting edge, (ii) a 1-vertex is an upper leaning point of some maximal 
segment but is not a 2-vertex, (iii) 0-vertices are all the remaining vertices. The 
number of i- vertices is denoted by n,. Given an orientation on the digital contour, 
the number of edges going from an z-vertex to a j-vertex is denoted by . 

Theorem 3. The number of supporting edges and of 1-vertices of F are related 
to its number of edges with 

"'^^^ <ni + 2n22. (5) 



i7(log to) 



An immediate corollary is that there are at least ne{F)/ f2 [log m) maximal seg- 
ments. 



Proof. From Theorem [2] and its Corollary [H we know that a DSS hence a max- 
imal segment cannot include more than n{logm) edges. Hence there cannot 
be more than f2(logm) 0-vertices for one 1-vertex or for one 2-vertex. We get 
'^oo < {ni + n2)^2(logm). We develop the number of edges with each pos- 
sible label: ne{r) = 7122 + ^02 + "12 + ?^20 + 't-21 + t^oo + t^oi + "-10 + 
Since, no2 + < 7122, + "21 < ^^22 and hqi + nio + rin < 3ni, we get 
ne{r) < 3n22 + "-00 + 3ni. Noting that a 2-vertex cannot be isolated by def- 
inition of supporting edges (Definition ^ gives 712 < 2ri,22. Once inserted in 
'^oo < ("-1 + n2)f^{logm) and compared with ne{r), we get the expected result. 

□ 

We now relate the DCP perimeter to the length of maximal segments. 

Theorem 4. The length of the smallest maximal segment of the DCP F is upper 
bounded: 

mmC\MSi) < n{\ogm)^^^. (6) 

I ne(r) 

Proof. We have Per{r) = J2n ^^{^i)- We now may expand the sum on sup- 
porting edges (22-edges), on edges touching a 1-vertex, and on others. Edges 
touching 1-vertices may be counted twice, therefore we divide by 2 their contri- 
bution to the total length. 



ni 



For the first term, each supporting edge indexed by j (a 22-edge) has an asso- 
ciated maximal segment, say indexed by j' . From Proposition [H we know that 
C\Ef)>\C\MSy). 

For the second term, each 1-vertex indexed by k is an upper leaning point 
of some maximal segment indexed by k' . Proposition [5] holds and C^[El^_^) + 
C\E]:^)>\C\MS,.). 

Putting everything together in Eq. |[7]), we get: 

Y,C\E,) > i^/:i(MS'y)+i^£i(Af5fcO > {^n22 + ^ni) mm C\M Si) > ^{n,+2n22) minC^MSi) 

"e >122 "1 

Inserting the lower bound of Theorem [3] into the last inequality concludes. 



4 Asymptotic properties of shapes digitized at increasing 
resolutions 

We may now turn to the main interest of the paper: studying the asymptotic 
properties of discrete geometric estimators on digitized shapes. We therefore 
consider a plane convex body S which is contained the square [0, 1] x [0, 1] 



(w.l.o.g.). Furthermore, we assume that its boundary 7 = 95 is with every- 
where strictly positive curvature. This assumption is not very restrictive since 
people are mostly interested in regular shapes. Furthermore, the results of this 
section remains valid if the shape can be divided into a finite number of convex 
and concave parts; each one is then treated separately. The digitization of S with 
step 1/m defines a digital convex polygon r{m) inscribed in a m x to grid. We 
first examine the asymptotic behavior of the maximal segments of r{m), both 
theoretically and experimentally. We then study the asymptotic convergence of 
a discrete curvature estimator. 

4.1 Asymptotic behavior of maximal segments 

The next theorem summarizes the asymptotic size of the smallest maximal seg- 
ment wrt the grid size to. 

Theorem 5. The length of the smallest maximal segment of r{m) has the fol- 
lowing asymptotic upper bound: 

mm {MS, {r{m))) < n{m^/^\ogm) (8) 

i 

Proof. Theorem[4]gives for the DCP r{m) the inequality min^ {M Si{r(m))) < 
i7(log to) ^n''{r{m)) ■ Since r{m) is convex included in the subset to x to of the dig- 
ital plane, its perimeter Per(_r(TO)) is upper bounded by 4to. On the other hand, 
Theorem [1] indicates that its number of edges ne{r{m)) is lower bounded by 
Ci(S')to^/'^. Putting everything together gives mini £^(M5'i(/^(TO))) < ^{^ogm) — 
which is once reduced what we wanted to show. □ 

Although there are points on a shape boundary around which maximal seg- 
ments grow as fast as 0{m^^^) (the critical points in [T2j), some of them do 
not grow as fast. A closer look at the proofs of Theorem |4] shows that a sig- 
nificant part of the maximal segments (at least J7(l/(logTO))) has an average 
length that grows no faster than f2{ni}^^ log to). This fact is confirmed with ex- 
periments. Fig. [21 left, plots the size of maximal segments for a disk digitized 
with increasing resolution. The average size is closer to m^^^ than to y/m. 

4.2 Asymptotic convergence of discrete geometric estimators 

A useful property that a discrete geometric estimator may have is to converge 
toward the geometric quantity of the continuous shape boundary when the dig- 
itization grid gets finer |3|4|10) . It may be expressed as follows. 

Definition 8. Let T he any geometric descriptor on the shape S with boundary 
7 and digitizations r{m). The discrete geometric estimator £ asymptotically 
converges toward the descriptor T for 7 iff 



\£{r{m)) - T{i)\ < e{m) with lim e{m) = 0. 

m— *-|-oc 



(9) 
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Fig. 2. For both curves, the digitized shape is a disk of radius 1 and the abscissa 
is the digitization resolution. Left: plot in log-space of the £^-size of maximal 
segments. Right: plot of the mean and standard deviation of the absolute error 
of curvature estimation, |k — 1| (expected curvature is 1). 



We now recall the definition of a discrete curvature estimator based on DSS 
recognition [3]. 

Definition 9. Let P be any point on a discrete contour, Q = B{P) and R = 
F{P) are the extremities of the longest DSS starting from P ('called half-tangents/ 
Then the curvature estimator by circumcircle is the inverse of the radius 

of the circle circumscribed to P, Q and R, rescaled by the resolution m. 

Experiments show that this estimator rather correctly estimates the curvature 
of discrete circles on average (w 10% error). It is indeed better than any other 
curvature estimators proposed in the litterature. Theorem B.4 of demon- 
strates the asymptotic convergence of this curvature estimator, subject to the 
hypothesis: 

Hypothesis 41 Half-tangents on digitized boundaries grow at a rate of f2{^/rn) 
with the resolution m. 

However, with our study of maximal segments, we can state that 

Claim. Hypothesis ST] is not verified for digitizations of C'^-curves with strictly 
positive curvature. We cannot conclude on the asymptotic convergence of the 
curvature estimator by circumcircle. 

Proof. It is enough to note that half-tangents, being DSS, are included in max- 
imal segments and may not be longer. Furthermore, since maximal segments 
cover the whole digital contour, some half-tangents will be included in the small- 
est maximal segments. Since the smallest maximal segments are no longer than 
fi(rn}/^ logm) (Theorem [5l, the length of some half-tangents has the same up- 
per bound, which is smaller than n[^Jm). □ 



The asymptotic convergence of a curvature estimator is thus still an open 
problem. Furthermore, precise experimental evaluation of this estimator indi- 
cates that it is most certainly not asymptotically convergent, although it is ac- 
tually on average one of the most stable discrete curvature estimator (see Fig. [21 
right). Former experimental evaluations of this estimator were averaging the 
curvature estimates on all contour points. The convergence of the average of all 
curvatures does not induce the convergence of the curvature at one point. 

5 Conclusion 

We show in this paper the relations between edges of convex hulls and maximal 
segments in terms of number and sizes. We provide an asymptotical analysis of 
the worst cases of both measures. A consequence of the study is the refutation 
of an hypothesis related to the asymptotic growth of maximal segments and 
which was essential in proving the convergence of a curvature estimator based 
on DSS and circumcircles [3]. Our work also applied to digital tangents since 
their convergence relies on the same hypothesis. The existence of a convergent 
discrete estimator of curvature based on DSS is thus still a challenging problem 
and we are currently investigating it. 
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A Proof of Proposition [2] 

Proposition [2l A pattern with an even complexity (say n = 2i) is such that 

UiLi = ((j2i-l+l,P2j-l-l) andLiU2 = (M2i-l)('72i-l,P2j-l) + (l?2i-2,P2i-2) + 

(—1,1). Moreover the DSS [Ui,Li\ has £^(^21-2)"^'"^ as a left factor, and the 
DSS [Li,U2] has i?(z2i_i)"^'~^ as a right factor. 

Proof. From Eq. ([3]) we have: P2i<l2i-i ~P2i-i'?2i = (—1)^*^^ = —1, which can be 
rewritten as: a(— 524-1) — b{—p2i-i) = 1 and eventually a{q2i — 92i-i) — b{P2i — 
P2i-i) = 1- We clearly obtain the Bezout coefficients. From its remainder we get 
the relatives coordinates of Li, as: UiLi = ((721-1 + l,P2j-i — 1)- From L1U2 = 
-UiLi + U1U2 we get : L1U2 = ((^21 - l)g2j-i + q2i~2 - 1, {_u2i - l)p2i-i + 
P2^^2 + 1). Using E{z2^) = Eiz2^-2)"''-'+^E{z2^-3)E{z2^-ir'-^-^ and UiLi = 
{u2i-iq2i-2 + q2i-z + l,U2t-iP2i-2 + P2i-3 " 1), it is clear that £;(z2i-2)"^'"' 
is a left factor of the DSS [J7iLi]. From Eq. |[2| and L1U2 we clearly see that 
£;(z2i_i)"^'-i is a right factor of the DSS [L1U2]. □ 

B Proof of Theorem H 

Lemma 4. We call P„ a pattern of complexity n whose Freeman code is E{zn). 
One can build strict right and left factors ( called respectively R and L) of Pn 
such that: 

(i) [RPn], [P-nL] and [RPnL] are DSS of slope z„, 
(ii) R and L are patterns (or successions of the same pattern) , 
(Hi) RPn, PnL and RPnL are not patterns, 

(iv) the slope of R is greater than that of Pn and the slope of Pn is greater than 



that of L, 

(v) maximal complexity of slope of R and L depends on parity of n: 



Complexity of Pn 


maximal complexity of R 


maximal complexity of L 


2i + l 


2i + l 


2i 


2i 


2i-l 


2i 



(vi) Complexity of factors obtained by substracting R or L from Pn depends on 



parity ofn: 



Complexity of P„ 


complexity of P„ \ R 


complexity of Pn \ L 


2i + l 


2i 


2i + 1 


2i 


2i 


2i-l 



Proof. Since R and L are strict factors of P„ , their Freeman moves are compat- 
ible with those of E{zn), giving same slope when i?,P„ and L are put together. 
Thus [RPnL] is a DSS of slope z„. From digital straightness we clearly have digi- 
tal convexity (see (5). Upper leaning points of this DSS are located at extremities 
of P„. 

We simply choose among strict factors R and L those that are patterns so 
that they fit descriptions given in Eq. ([T| and Eq. We may now describe 
them given the parity of n. 

Consider the case where n is odd (say n — 2i + 1), from Eq. H)) we get: 
R = E{z2i)'''''+'-''E{z2i-i) and L = £;(z2»)"''+'"' with r > and Z > 0. If i? and 
L are longer patterns, they are not anymore strict factors of P2i+i • We see that R 
is a pattern of complexity 2i + l and that L is a succession of the pattern E{z2i), 
with a complexity of 2i. The slope of i? equals ^22+1 ~ [0, ui, . . . , U2i,U2i+i—r] = 

g±i. FVom Eq. Q we get that ^ = The sign of z',^^, - is 

that of P2j_|_ig2i — q'2i+iP2i, and is positive (see Eq. ([3])). Thus the slope of R is 
greater than that of P2i+i. Same reasoning applied to Z2i+i — Z2i brings that the 
slope of P2i-i-i is greater than that of L. Factor obtained by substracting R from 
P2i+i equals E{z2iY and substracting L from P2i+i gives i?(z2i)'^(^2i-i)- 

Consider now that n is even (say n = 2i), from Eq. ^ we get: R = 
E{z2i-i)^^'~^ and L = E{z2i-2)E{z2i-iY'^'~^ ■ If R and L are longer patterns, 
they are not anymore strict factors of P2i ■ Clearly, R has a complexity of 2i — 1 
and that of L equals 2i. The slope of L equals zL = [0, wi, . . . , U2i-i,U2i—l] = 

From Eq. (g]) we get that ^ The sign of Z2i - 4, is that of 

q'2iP2i-i — p'2il2i-i, and is positive (see Eq. Thus the slope of P„ is greater 
than that of L. Same reasoning applied to Z2i-i — Z2i brings that the slope of R 
is greater than that of P„. Factor obtained by substracting R from P2i equals 
E{z2i-2)E{z2i-iY and substracting L from P2i gives E{z2i-iy ■ 

It is now clear that slopes are strictly decreasing from R to P„ and from P„ 
to L whatever the parity of n. □ 

Theorem 6. Let E be a supporting edge whose slope has a complexity n, n >2, 
then the maximal segment containing E includes at most n other edges on each 
side of E. 

Proof. We construct 2n digital edges around E: 

- iRi)i<i<n at left of E, 

- {Li)i<i<n at right of E. 

These edges are such that [P„ . . . Ri . . . RiELi . . . Lj . . . P„] is a DSS of slope 
Zn = a/b and has no other upper leaning points but those located on E. E may 



contain several times the pattern E{z„). It is clear that Rn ■ ■ ■ Ri ■ ■ ■ Ri (resp. 
Li. ..Lj . . .Ln) has to be a right (resp. left) strict factor of E{zn). Moreover 
Ri is a right strict factor of E{zn) \ i?i . . . Ri-i and Li is a left strict factor of 
E{zn) \ ii . . . Li-i. From Proposition [3] any of the digital edges {Ri)i<i<n and 
(Li)i<i<„ is a pattern or a succession of the same pattern. From Eq. |(T|) and 
Eq. ^ two successive digital edges with same complexity (say n) cannot form a 
right or left strict factor of a pattern with same complexity. Thus complexities 
of {Ri)i<i<n and (ii)i<i<n are decreasing when i increases. Moreover to fullfil 
convexity properties, slopes of edges are decreasing from i?„ to L„. 

We now build {Ri)i<i<n when n is odd (say n = 2i + 1). From Lemma[H Ri 
has a complexity that equals 2i + 1 and i?2 is a right strict factor of a pattern 
whose complexity equals 2i. Applying Lemma H] brings R2 with a complexity 
of 2i — 1 . Applying the same reasoning recursively brings other edges as shown 
on Table [TJ Lemma |4] also implies decreasing slopes and give upper bounds in 
complexity of factors. 

Constructions for the three other cases are given in Tables [T] and [2] and 
follow the same reasoning. To satisfy full decomposition each {uk)i<n has to 
be equal or greater than 2. If this condition is not meet for some k, than steps 
associated with it (e.g. any factors containing Uk — Vj or Uk — lj as powers of 
some pattern) are skipped. This concludes the proof. □ 



Table 1. Constructions of {Ri)i<i<n and (ii)i<i<„ given n odd . 



Constructions of (-Ri)i<i<n when n = 2i + 1 



Factor 


Freeman moves 


Complexity 


Ri 


E{z2^r^'+'-'^^E(z2^-l) 


2i + 1 


R2 


E{Z2^^^Y^^-''^ 


2i - 1 


Rs 




2i - 1 


Ri 


E{Z2,-3Y^^-^-'-^ 


2z - 3 








R2j 


E{z2,+ l-2jr^^+^-^'-''^' 


2i + 1 - 2j 


R-2i + l 


E{z2^-2jY^^+^-^' + ' E{z2i-l-2j) 


2i + l- 2j 








R2i+1 


QUl-rai+i ^ 


1 



Constructions of (I/i)i<i<Ti when n = 2i + 1 



Factor 


Freeman moves 


Complexity 


Li 




2i 


L2 


E{z2^-^2)E{z2^^lr^'-'■' 


2i 


L3 


E{z2i-2r'^-'-" 


2i-2 


U 


E(z2^-i)E{z2r-3T''^-^-'^ 


2i-2 








L2] 


E{z2^-2J)E{z2^ + 1^2JT^' + ^-^=''^-' 


2i + 2- 2j 


L2j + l 




2i - 2j 








L2i+1 


QUI — i2i+l 






Table 2. Constructions of {Ri)i<i<n and (ii)i<i<n given n even. 



Constructions of {Ri)i<i<n when n = 2i 



Factor 


Freeman moves 


Complexity 


Ri 




2i - 1 


R2 




2i-l 


Rs 




2i - 3 


Ri 


E{z2^-4T^^-'-''^E{z2^-5) 


2i - 3 








R2j 


E{z2i-2jT^'+'-^'-''^'E{z2,-l-2j) 


2i + 1 - 2j 


R-ij+i 


E{z2i-l-2jT^'-^=-''^= + ^ 


2i - 1 - 2j 








R2, 


QUl-r2i ^ 


1 



Constructions of (I/i)i<!<n when n — 2i 



Factor 


Freeman moves 


Complexity 


Li 


E{z2r-2)E{z2^-lr^'^'' 


2i 


1.2 


£;(22«-2)"^'-l-'= 


2i-2 


La 


E{z2^-4)E{z2^-3r^'-'''^ 


2i-2 




E{z2i-4r^'-^-'^ 


2i - 4 








L2j 


E{z2r-2jY^^ + ^-^'-'^' 


2i - 2j 


L2j+1 


i;(z2,-2-2j)S(z2«-l-2j)"^'-=^-'=^ + l 


2i - 2j 








L2i 








